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Abstract. The Narrow Line Seyfert 1 galaxy Mrk 766 
has been observed by BeppoSAX on May 1997. The source 
was fairly variable, both in flux and spectrum, during the 
observation. The variability is the largest around 2 keV, 
possibly due to variations in the warm absorber proper- 
ties. An oxygen line is clearly detected, at least in the first 
part of the observation. An edge around ~7.5 keV, coupled 
with the lack of any detectable iron line, suggests either 
reprocessing of the primary X-rays by a mildly ionized 
disc or absorption by a further, thicker and more ionized 
material. 

Key words: Galaxies: individual: Mrk 766 - Galaxies: 
Seyfert - X-rays: general 



1. Introduction 

Mrk 766 is a rather bright (typical flux of about one milli- 
Crab) Narrow Line Seyfert 1 galaxy at redshift 0.012, 
showing complex and puzzling X-ray variability (Molendi, 
Maccacaro & Schaeidt 1993; Molendi & Maccacaro 1994; 
Page et al. 1999) on time-scales as short as ~1000 sec 
(Leighly et al. 1996). The 1992 June ROSAT PSPC ob- 
servations, discussed by Molendi & Maccacaro (1994), in- 
dicated large amplitude variations in the softest part of 
the spectrum (0.1-0.9 keV), while the hardest part (0.9-2 
keV) was consistent with being constant. The authors sug- 
gested accretion disc emission as the origin of the variable 
component. Recently, Page et al. (1999) analysed spectral 
variabihty in ROSAT PSPC data collected from 1991 to 
1994, reaching different conclusions: a highly variable (on 
time-scales as short as 5000 s) hard power law component, 
and a constant soft excess. The ASCA-ROSAT observa- 
tion reported by Leighly et al. (1996) shows a complex 
and variable spectrum, with the primary power law pho- 
ton index changing from 1.6 to 2 within the observation, 
and variable warm absorber too. The soft excess variabil- 
ity was anti-correlated with that of the hard component. 



A narrow iron Ka line was detected only when the source 
was in a high flux level. 

We observed Mrk 766 with BeppoSAX to take advan- 
tage of the wide energy band of the satellite and try to 
better understand the origin of the spectral complexity in 
this source. 

The paper is organized as follows: Sec. 2 describes the 
observation and data reduction; Sec. 3 presents the tem- 
poral, and Sec. 4 the spectral, analysis; the results are com- 
pared with the ASCA-ROSAT ones in Sec. 5, discussed 
in Sec. 6 and summarized in Sec. 7. 



2. Observation and data reduction 

BeppoSAX (Boella et al. 1997) observed Mrk 766 on 17-18 
May, 1997 as part of a Core Program devoted to a spectral 
survey of bright Seyfert 1 galaxies. The exposure time was 
79 ksec for the MECS, the HPGSPC and the PDS and 42 
ksec for the LECS. We discuss here data from the LECS, 
MECS and PDS only, because the source is too faint for 
the HPGSPC to be profitably used. 

Spectra and light curves were extracted from the imag- 
ing instruments within circles of 8' (LECS) and 4' (MECS) 
radii centred on the source. Background spectra taken in 
the same regions from blank sky observations have been 
used for subtraction. Regarding the PDS, background sub- 
tracted spectra automatically generated at the SAX Scien- 
tific Data Centre (SDC) have been used. Fixed Rise Time 
thresholds have been adopted. 

Time-averaged count rates are: 0.2765±0.0026 (LECS, 
0.1-10 keV); 0.2624±0.019 (MECS, 1.5-10 keV, 2 unitfQ); 
0.2863±0.0038 (PDS, 15-200 keV). The 2-10 keV mean 
flux (obtained with model 5 in Tables 1 and 2) is 
2.05x10"" ergcm-2 g-i. 

In the field of view of the PDS there are two possible 
confusing sources, both of them BL Lac objects. At an 
angular distance of about 20 arcminutes there is ON 325. 



^ The observation was performed after the failure of one of 
the three MECS units. 
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The source is clearly visible in the MECS, with a count 
rate of about 1/20 that of Mrk 766. The faintness of the 
source and its off-axis position do not permit a detailed 
spectral analysis, but the MECS spectrum appears to be 
rather steep, and therefore its contribution to the PDS 
is likely to be very small. More serious is the case of the 
other source, 2A 1219+305 (PG 1218+304), a member of 
the Piccinotti sample (Piccinotti et al. 1982) with a mean 
2-10 keV flux of 3x10"" erg cm^^ s"^ and a 2.5 a de- 
tection with BATSE (Mahzia et al. 1999). The source is 
about 44' arcminutes away from Mrk 766, i.e. outside the 
MECS field of view, and it is therefore not possible to 
estimate its contribution to the PDS, which may be sub- 
stantial. For this reason we preferred to exclude the PDS 
from the spectral fitting, using the data a posteriori only 
as an upper limit to the hard X-ray emission of Mrk 766. 

3. Data analysis and results. Temporal analysis 

The LECS (0.1-10 keV) and MECS (1.5-10 keV) light 
curves are shown in Fig. n^. The source varies as much 
as a factor ~2 on time-scales of a few thousands of sec- 
onds. To search for spectral variability, we have also plot- 
ted the (1.5-4)7(0.1-1.5) LECS hardness ratio, and the 
(4-10)7(1.5-4) MECS hardness ratio. Clearly, there is sig- 
nificant spectral variability, especially in the second part 
of the observation. To investigate the energy dependence 
of this variability, we calculated the Normalized Excess 
Variance (NEV; see Nandra et al. 1997) (Fig. |) for differ- 
ent energy intervals, using a 128 s binning. The source is 
variable at all energies, but the amplitude is the highest 
around 2 keV. 

The PDS light curve and the PDS/MECS ratio are 
also shown in Fig. |l|. The light curve is very variable, the 
background subtracted PDS count rates sometimes going 
down to values consistent with zero. The light curve is 
different from those of the LECS and MECS, suggesting 
either a strong spectral variability or a significant contam- 
ination by the nearby BL Lac, 2A 1219+305. 

4. Data analysis and results. Spectral analysis 

The spectral analysis was performed on the LECS and 
MECS instruments only, because of the possible confusion 
problems in the PDS outlined above. The PDS data have 
been used only to check a posteriori that a LECS+MECS 
best fit spectrum, when extrapolated to higher energies, 
would not exceed the observed data. 

All fits described below have been performed with the 
XSPEC software package. Quoted errors refer to 90% confi- 
dent level for two interesting parameters (i.e. Ax^=4.61). 

Because of the significant spectral variability, it is not 
wise to use the spectrum averaged over the whole obser- 
vation. On the other hand, it is important to collect as 
many photons as possible to search for spectral details. 
As a trade-off we have divided the observation in two 
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Fig. 1. From top to bottom: LECS light curve and hard- 
ness ratio; MECS light curve and hardness ratio; PDS light 
curve (background subtracted) and PDS/MECS hardness 
ratio. 



parts, and extracted spectra from: a) the beginning of the 
observation till 8x10* sec (see Fig. |l|), because within this 
interval variations in the hardness ratio, especially in the 
MECS, are not dramatic; b) from this time till the end of 
the observation. 

As it will be seen in the following (Tables 1 and 2), no 
models will give a fully acceptable fit, the null hypothesis 
probability being 0.12 at most. This is likely due to the 
spectral variability discussed in the previous section, as 
our time selection makes this problem alleviated but not 
completely cured. 

4-.1. First half of the observation 

Let us start discussing the first half of the observation. 
A simple power law absorbed by the Galactic column (i.e. 
1.77x10^° cm"2; Elvis, Wilkes & Lockman 1989) is a very 
poor fit to the data (see Fig. and Table 1, model 1), 
with an unacceptable reduced x^ of 1.70. From the fig- 
ure, it is clear that most of the contribution to the x^ 
comes from the softest part of the spectrum, i.e. below ~1 
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Fig. 2. The Normalized Excess Variance for different en- 
ergy intervals. 



keV, while no evidence for typical Seyfert 1 components 
like the iron line or the reflection continuum is present. A 
deficit of counts above 7 keV is also apparent. Leaving the 
column density free to vary (model 2), the quality of the 
fit significantly improves, but the fit is still unacceptable 
(X?=1.50). 

4.1.1. The warm absorber 

Inspection of the residuals suggests the presence of warm 
absorption features. We therefore added the ABSORI model 
(model 3 in Table 1), obtaining a significant improvement 
in the quality of the fit. The warm absorber material has 
been assumed to be ionized by the observed power law. 
The best fit value of the ionization parameter, ^, is about 
240 erg cm s^^; with this value, the most prominent edge is 
the O VIII one. The column density of the warm absorber 
is ~8xl02i cm-2. 



4.1.2. The Oxygen line 

A further inspection of the residuals reveals an excess 
around 0.6 keV, which can be well fitted (model 4) by 
a narrow gaussian line at 0.59(±0.05) keV (source rest 
frame), with an equivalent width of 68(l4g) cV. The sta- 
tistical significance of this line is 97.3% (F-test). The line 
energy is consistent with He-like oxygen (0.57 keV); a fit 
with the line fixed at this energy, plus a line at 0.65 keV 
(H-like oxygen) does not provide a statistical improve- 
ment; the best fit EWs are 60 and 30 eV, respectively. 



Fig. 3. The spectrum along with the best fit model (upper 
panel) and data/model ratio (lower panel) for the first half 
of the observation, when fitted with a simple absorbed 
power law. 



The value of ^ is now lower. The best fit ionization struc- 
tures of the absorber and the emitter are somewhat dif- 
ferent (the most relevant ion being O vii in the emitter, 
and O VIII in the absorber), but still consistent each other 
within the (fairly large) errors. Finally, a relativistic line 
(diskline model in xspec) fits the line equally well as a 
narrow gaussian line. Fixing the inner radius at 6 r^, i.e. 
the last stable orbit for a Schwarzschild black hole, we ob- 
tain an outer radius of thousands of Vg and an inclination 
angle consistent with zero (i.e. face-on disc). 



4.1.3. The iron edge. Ionized reflection of ionized 
absorption? 

While the above fit is satisfactory from a statistical point 
of view, a deficit of counts above 7 keV is still apparent. 
Adding an absorption edge to model 4 the fit actually 
improves, the best fitting parameters being an edge en- 
ergy of 7.55 keV, and t=0.26. This edge may be related 
to reflection from circumnuclear matter, and therefore we 
added a Compton reflection component. We allowed the 
matter to be ionized (pexriv model in xspec), both be- 
cause the edge is at energies larger than 7.1 keV, the value 
of neutral iron, and because the lack of any observed iron 
Ka line (see next paragraph) is only possible, in presence 
of a Compton reflection component, if the iron is mildly 
ionized and therefore resonant destruction possible (e.g. 
Matt, Fabian & Ross 1993, 1996). The results are reported 
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in Table 1 as model 5. Interestingly, the value of the ion- 
ization parameter of the reflecting matter, even if poorly 
determined, is consistent with what is expected in order 
to have the iron line significantly destroyed. It is worth 
noticing that such an ionized disc could also account for 
the observed oxygen emission line, even if the equivalent 
width is larger than expected by a factor of a few, pointing 
to a possible oxygen overabundance. Moreover, the best 
fit parameters of the warm absorber material are now such 
that the most important oxygen ion is the Helium-like. It 
is therefore possible that both the accretion disc and the 
warm absorber contribute to the line flux. 

Alternatively, the iron edge may be due to absorbing 
material. We therefore added a second ionized absorber 
(model 6), instead of the reflector. The fit is as good as 
the one with the reflector. This second absorber results 
to be more thick and ionized than the other one (but the 
parameters are loosely constrained) . It may be interesting 
to note that a similar double-absorber solution has been 
found in the (broad line) Seyfert 1 NGC 3516 (Costantini 
et al. 2000). 

4.1.4. The iron line 

Leighly et al. (1996) detected a narrow 6.4 keV iron line 
with an equivalent width of about 100 eV in the ASCA 
spectrum of Mrk 766, but only when the source was in a 
high state. We searched for an iron line in the BeppoSAX 
data, but could find only upper limits. The upper limits to 
a narrow iron line are 110, 43 and 36 eV if the line energy 
is fixed to 6.4 keV (neutral iron), 6.7 keV (H-like iron) or 
6.97 keV (H-like iron), respectively. Therefore, our result 
is marginally consistent with the line found by Leighly et 
al. (1996). Of course, the upper limits would increase if 
the line is broad. As explained in the previous section, 
the lack of an observable line is still consistent with the 
reflection scenario, because the best fit ionization state is 
such to have the line destroyed by resonant re-absorption. 

4.1.5. PDS data 

As shown in Fig. 0, the extrapolation of model 5 to the 
PDS data is rather good, which may be an indication (but 
unfortunately not a proof) that the contribution to the 
PDS count rate from 2A 1219+305 is low. Also the ex- 
trapolation of model 6 does not exceed the observed PDS 
data. 

4.1.6. The soft excess 

We then searched for a soft excess (modeled as a black 
body), a component which has been already observed 
in this source (see Sec. 1), and whose presence may be 
expected by analogy with other Narrow Line Seyfert 1 
Galaxies (e.g. Boiler, Brandt & Fink 1996). The fit with 
the soft excess instead of the warm absorber (other com- 




1 10 

Energy (keV) 



Fig. 4. The spectrum and best fit model (upper panel) 
and residuals (lower panel) for model 5 in Table 1. The fit 
has been performed on the LEGS and MEGS data only, 
and then the best fit model has been extrapolated to the 
PDS data. 

ponents as in model 5) gives a worse Xr (1-23/119; best 
fit temperature of about 30 eV), while the addition of the 
black body to model 5 does not provide a significant im- 
provement (Xr=l-13/117; best fit temperature of 200 eV). 
We therefore conclude that such a component is not re- 
quired by the data, and confirm that the warm absorber 
is real, and not an artifact of not having included the 
soft excess. However, it should be noted that while the 
fits presented in Table 1 have been obtained with a col- 
umn density of the cold absorber in excess of the Galactic 
value, the column density turns out to be close to the 
Galactic one when the black body component is included. 
Therefore, the presence of the soft excess remains rather 
ambiguous. It is worth noting that the BeppoSAX power 
law index is significantly steeper than the one observed 
by ASGA, while the average 2-10 keV flux is about the 
same in the two observations, making any soft excess in 
the BeppoSAX data difficult to observe. 

4-2. Second half of the observation 

We then analysed the spectra extracted from the second 
half of the observation. For the sake of conciseness, only 
models 5 and 6 are reported in Table 2. The oxygen line 
is now not required by the data, and we fixed its energy 
to 0.6 keV (see Table 1; this energy would correspond to a 
blend of He- and H-like atoms) to get an upper limit. The 
fit with two absorbers is now preferable to that with one 



Table 1. Best fit parameters for the first half of the observation 
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" Column density (in units of 10 cm ) of the second absorber in model 6 



absorber and one reflector, from a statistical point of view. 
Again, a blackbody component instead of the warm ab- 
sorber gives a much worse fit. This time, however, adding 
the blackbody component a significantly better fit is ob- 
tained, but at the expense of a very steep power law com- 
ponent (r=3.4) and an unplausibly large refiection com- 
ponent {r~21). As there is significant spectral variations 
during the second half of the observation, it is possible 
that this improvement occurs because a further compo- 
nent may help fitting the fictitious spectral complexity 
arising from time-averaging over different spectral states. 
Not surprisingly, the x^ is significantly higher than in the 
first half of the observation. 

Let us now compare the results obtained in the two 
halves of the observation. Two main differences are evi- 
dent: first, in the second half of the observation the oxygen 
line is not required by the data, and only an upper limit on 
its EW can be obtained. This limit, however, is consistent 
with the value for the first half. Second, the WA (both of 
them in model 6) appears to be thicker and more ionized 
in the second half than in the first half (even if, within the 
errors, the values of the ionization parameters are consis- 
tent with each other) . The two best fit absorption models 
are shown in Fig. g[ where it can be seen that the dif- 
ferences between the two models are large in the ~0.9-2 
keV interval. Intriguingly, these are the energies where the 
NEV has a maximum (Fig. ^). Unfortunately, the limited 
statistics prevent from a more detailed analysis. 



5. Comparison with previous observations 

The behaviour of Mrk 766 during the BeppoSAX obser- 
vation is rather different than during the ASCA-ROSAT 
observation, performed on December 1993 and discussed 




Energy 



Fig. 5. Best fit warm absorber model for the two halves of 
observations, obtained with model 5 in the Tables. Solid 
line: first half. Dash-dotted line: second half. 



by Leighly et al. (1997). The source was at an average 
flux level very similar in the two observations. While in 
both cases there is significant fiux and spectral variability 
on time scales of thousands of seconds, the origin of the 
spectral variability appears to be different. In the ASCA- 
ROSAT observation, it is driven by significant changes in 
the spectral index of the power law, while in our obser- 
vation it seems to be mostly related to changes in the 
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Table 2. Best fit parameters for the second half of the observation 
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warm absorber properties. Indeed, in our observation the 
power law index is consistent with being constant, but at a 
value significantly steeper (i.e. T ~2.2) than those found 
by Leighly et al. (F ~1.6-2). Moreover, in the ASCA- 
ROSAT observations a soft excess is strongly required by 
the data, while in our observation this component may 
be absent. (It is worth noting than even in the ASCA- 
ROSAT observations the soft excess is marginally visible 
when the power law is the steepest, because the soft com- 
ponent is almost constant during the observation; it may 
therefore be not surprising that in our observation the 
soft component is hardly detectable). A further difference 
is the presence, in the ASCA data and during the high 
flux phase, of an iron line with ^100 eV equivalent width, 
while in the BeppoSAX data such a line is not visible (but 
the upper limit on the equivalent width is marginally con- 
sistent with the abovementioned value). 

To summarize, in both observations the source under- 
went large and complex variability, but the pattern and 
origin of the variations seem to be very different in the 
two cases. 

6. Discussion 

While the Narrow Line Seyfert 1 Galaxy Mrk 766 shares 
with other members of its class the large amplitude and 
short time-scale variability, it does not share (at least dur- 
ing the observation described here) another characteristic 
which is common in Narrow Line Seyfert 1 galaxies, i.e. a 
prominent soft X-ray emission. The power law component 
is, actually, somewhat steeper than the typical value for 
classical Seyfert 1 galaxies, but the soft excess is modest, 
if present at all. The warm absorber is also typical of clas- 
sical, 'broad lines' Seyfert Is, being dominated by oxygen 
and neon edges. No ~1 keV feature (besides those related 
to the warm absorber and cured by the inclusion of this 
component) is present, differently from what observed in 
other sources of this class (e.g. Leighly et al. 1997; Fiore 
et al. 1998; Leighly 1999; Vaughan et al. 1999a). This is 
not surprising, as these features have been observed so far 
only in very steep spectrum sources, a fact naturally ex- 
plained if these features are due to a blend of resonant ab- 



sorption lines, mainly from iron L shell (Nicastro, Fiore & 
Matt 1999). The column density and ionized parameter of 
the warm absorber appear to have been changed between 
the two halves of the observations, in agreement with the 
fact that the variability is the largest around 1 keV. A so 
dramatic change of the column density, however, is un- 
plausible, and it may be an artifact of having used pure 
photoionization equilibrium, single-zone models, while in 
reality it is possible that: the absorber is out of equilibrium 
most of the time (not surprisingly, given the flux variabil- 
ity; see Nicastro et al. 1999 for non equilibrium models); 
there is some contribution from coUisional ionization; the 
absorbing medium is geometrically or physically complex. 

Similar to other Narrow Line Seyfert Is (TON S 180: 
Comastri et al. 1998, Turner, George & Nandra 1998; Ark 
564: Vaughan et al. 1999b) is, instead, the possible pres- 
ence of reprocessing from ionized matter. Ionized accretion 
discs are expected when m is high (Ross & Fabian 1992; 
Matt, Fabian & Ross 1993), and accretion rates close to 
the Eddington one have been indeed invoked to explain the 
NLSl phenomenon (e.g. Pounds, Done & Osborne 1995; ). 
The (mild) ionization of the reprecessor may also account 
for the lack of any observed iron line (Matt et al. 1993; 
1996), as the line photons may be resonantly trapped and 
then destroyed by the Auger effect. 

Reflection from a mildly ionized disc may also explain 
at least part of the O Vil emission line observed in the first 
half of the observation (in the second half only an upper 
limit is obtained, which however is consistent, within the 
error, with the value measured in the first half). In this 
case, the line is expected to be broadened by relativistic 
effects; the quality of our data is not good enough to dis- 
tinguish between a narrow and a relativistic line, provided 
that the inclination angle of the disc is low. An oxygen line 
possibly from an ionized disc was observed by Piro et al. 
(1997) in the ASCA data of E1615-h061, while the same 
line observed by BeppoSAX in the spectrum of NGC 5548 
(Nicastro et al. 2000) clearly originates from outflowing 
material, as demonstrated by the Chandra/LETG obser- 
vation (Kaastra et al. 2000). In Mrk766, it is possible that 
both the accretion disc and the warm absorber contribute 



to the observed emission, as the best fit ionization struc- 
ture for both materials includes a significant fraction of 
O VII. High energy resolution observations are needed to 
definitely settle this issue. 

7. Summary 

The main results of the BeppoS AX observation of Mrk 766 
may be summarized as follows: 

• The spectrum of Mrk 766 is well fitted by an ab- 
sorbed power law (with cold absorption in excess of the 
Galactic one), a warm absorber, an oxygen line and a 
further component (mainly accounting for the iron edge) 
which may be either a second and more ionized warm ab- 
sorber or a mildly ionized disc. A soft excess is not required 
by the data, but its inclusion makes the cold absorber 
closer to the Galactic value. 

• The spectral variability during the observation may 
be entirely due to variations in the properties of the warm 
absorbcr(s), even if the behaviour of the absorber is com- 
plex and no obvious correlation with the flux is apparent. 

• The power law is significantly steeper than during 
the ASCA-ROSAT observation (Leighly et al. 1996). The 
spectral variability pattern is also very different, indicat- 
ing that we have observed the source in a different state. 

• There is clear evidence for an O vii emission line, 
which may originate cither in the accretion disc or in the 
warm absorber, or, most plausibly given the large equiva- 
lent width, in both. 
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